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Abstract 

We establish by means of a computer search that a complete graph 
on 14 vertices has 98,758,655,816,833,727,741,338,583,040 distinct and 
1,132,835,421,602,062,347 nonisomorphic one- factorizations. The enu- 
meration is constructive for the 10,305,262,573 isomorphism classes 
that admit a nontrivial automorphism. 
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1 Introduction 

A one-factor of a graph is a 1-regular spanning subgraph. A one-factorization 
of a graph is a set of one-factors such that every edge of the graph occurs 
in a unique one-factor. Two one-factorizations are isomorphic if there is a 
bijection between the vertex sets of the underlying graphs that maps the 
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one-factors in one one-factorization onto the one-factors of the other. One- 
factorizations of complete graphs are basic combinatorial objects with a wide 
variety of applications [U [HI [19]. 

A fundamental problem for any family of combinatorial objects is that of 
classifying these for small parameters, that is, of listing exactly one represen- 
tative from each isomorphism class. The complete graph K n on n vertices 
has, up to isomorphism, a unique one-factorization for n = 2, 4, 6. The graph 
K$ has 6 nonisomorphic one-factorizations; this result was obtained by Dick- 
son and Safford [2] almost exactly one century ago. In 1974, Gelling and Odeh 
[I] published the result that K\q has 396 nonisomorphic one-factorizations. A 
result by Dinitz, Garnick, and McKay [3] showing that there are 526,915,620 
nonisomorphic one-factorizations of Ki 2 appeared exactly two decades af- 
ter the result by Gelling and Odeh. The result for Kyi was one of the 
most extensive computer classifications at that time, and — with an estimated 
1.132 x 10 18 nonisomorphic one-factorizations [3 J of K14 — the current limit 
for methods that rely on constructively listing one specimen from each iso- 
morphism class was thereby reached. 

For many types of combinatorial objects, no essentially faster method 
than explicitly constructing one representative from each isomorphism class 
is known for counting the number of isomorphism classes. However, in cases 
where one is able to count the "labeled" objects substantially faster, the 
following alternative approach for determining the number of isomorphism 
classes becomes possible. 

Let T be a finite group that acts on a finite set Q. Denote by Ni the 
number of orbits on Q whose elements have stabilizer subgroups of order i in 
T. Then, by the Orbit-Stabilizer Theorem, 

mi = |r|£f- (i) 

i 

If we know |T|, |0|, and N for each % > 2, then we can solve (JTJ for N\ and 
thereby obtain the number of orbits ^ N. This observation can be applied 
in a combinatorial context by considering a group action for which (i) the 
set Q consists of the "labeled" objects; (ii) the orbits of T on Q correspond 
to the isomorphism classes of objects (the "unlabeled" objects); and (iii) the 
stabilizer subgroup of an object under the action of T corresponds to the 
automorphism group of the object. In particular, the values Ni for i > 2 are 
obtained by classifying up to isomorphism the objects that have a nontrivial 
automorphism group. 
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Latin squares and one-factorizations of complete graphs are two types of 
objects for which the outlined idea is effective. Recently, McKay, Meynert, 
and Myrvold [H] successfully applied this general idea to the problem of 
enumerating Latin squares. In the present paper, it is applied to the problem 
of enumerating one-factorizations of the complete graph K14. 

The approach for counting the number of distinct one-factorizations of 
K u is discussed in Section^ there are 98,758,655,816,833,727,741,338,583,040 
such objects. Classification of one-factorizations of Ku with a nontrivial 
automorphism group is considered in Section [31 Application of p]) then 
reveals that there are 1,132,835,421,602,062,347 isomorphism classes of one- 
factorizations of K14. Particular emphasis is put on verifying the correctness 
of these computational results. 

Conventions. For standard graph-theoretic terminology we refer to [20J. 
All graphs considered in this paper are undirected and without multiple edges 
or loops. For a graph G, denote by G the complement of G, by T{G) the set 
of all one-factors of G, by [G] the isomorphism class of G, by Aut(G) the au- 
tomorphism group of G, by LF(G) the number of distinct one-factorizations 
of G, and by NF(G) the number of nonisomorphic one-factorizations of G. 

2 Counting distinct one-factorizations 

To count the number of distinct one-factorizations of K u , that is, to compute 
LF(ifi4), we essentially rely on a recursion pioneered by Dinitz, Garnick, 
and McKay (3j Sect. 4], who used the recursion to check their constructive 
enumeration of the isomorphism classes of one-factorizations of Kyi- 

In our case, however, we do not have an independently computed value of 
LF(i^ 14 ) available. To gain confidence that the computed value of LF(iT 14 ) is 
correct, we will use (i) a new "forward accumulation" technique to compute 
the value of LF(i£i4); and (ii) Dinitz-Garnick-McKay recursion to check 
both the intermediate results and the final computed value. 

It is convenient to start by describing Dinitz-Garnick-McKay recursion 
and then proceed to describe the new technique. 

2.1 Dinitz— Garnick— McKay recursion 

Consider a graph G that has a one-factorization. Such a graph is necessarily 
/c-regular for some nonnegative integer k. Avoiding trivial cases, we assume 
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that k > 2. Let F be a one-factor in G, and denote by G — F the (A; — 1)- 
regular graph obtained by deleting the edges of F from G. To arrive at 
a recursion, suppose that we know the value of LF(G — F) for every one- 
factor F of G. By counting in two different ways the number of distinct 
one-factorizations of G with one individualized one-factor, we have 

k ■ LF(G) — LF(G-F). (2) 

F&F{G) 

In particular, if we know LF(if) for every (k — l)-regular graph H, then 
we can compute LF(G) for any given fc-regular graph G via ([2]). This is 
Dinitz-Garnick-McKay recursion from {k — l)-regular to fc-regular graphs. 

In practice it suffices to evaluate the right-hand side of (J2J) only for one- 
factors that contain a fixed edge of G, in which case the multiplication by k 
on the left-hand side is not necessary. Similarly, LF(if) needs to be computed 
for only one graph H in each isomorphism class of regular graphs. 

Dinitz-Garnick-McKay recursion has the property that it considers each 
k- regular graph in turn, and "looks back" at the (k — l)-regular graphs. We 
describe next a new technique that "looks forward" at the /c-regular graphs 
from each (k — l)-regular graph in turn. 

2.2 Forward accumulation 

The forward accumulation approach is based on the immediate observation 
that every fc-regular one-factorizable graph G can be decomposed into a 
union G = H U F of (i) a (k — l)-regular one-factorizable graph H; and (ii) 
a one-factor F G J-(H). Put otherwise, we can visit every isomorphism class 
of /c-regular one-factorizable graphs by the following procedure: for each 
isomorphism class [H] of (k — l)-regular one-factorizable graphs, consider 
exactly one graph H from the isomorphism class; for each such graph H, 
consider each one- factor F G J-(H); for each such pair (H,F), visit the 
isomorphism class [H U F] . 

To compute the value LF(G) for each visited isomorphism class [G], we 
associate with [G] an accumulator variable xug\ that is initially set to zero 
and incremented whenever [G] is visited. Our objective is to have the value 
k-LF(G) in the accumulator when the visiting procedure halts. To determine 
an appropriate increment to x\g\ on each visit, we proceed to analyze the 
visiting procedure in more detail. 
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To this end, consider the set of all pairs (H, F) such that H is a (k — Ir- 
regular one-factorizable graph and F G F{H). Let us view two such pairs as 
isomorphic if one can be obtained from the other by relabeling the vertices. 
The following lemmata are immediate consequences of the Orbit-Stabilizer 
Theorem. 

Lemma 1 Any graph G in the class [H U F] admits exactly 

a(H F)- < Aut ( gUF )l 
1 ' ; |Aut(#)nAut(F)| 

decompositions G = H' U F' into pairs (H', F') in the class [(H, F)} . 
Lemma 2 The procedure visits a class [H U F] exactly 

y ' ; |Aut(£T)nAut(F)| 
times via pairs (H',F') in the class [(H,F)]. 

It now follows from Lemma [1] and ([2]) that, for any /c-regular graph G, 
v(H,F).LF(H)= J2 LF(H) = fc.LF(G). (3) 

[(H,F)]:HUF=G (H,F):HUF=G 

This observation enables us to accumulate the value k ■ LF(G) for each k- 
regular G. Namely, each time [G] is visited via a pair (H,F), we increment 
X[a\ by the rule 

x [G ] <- X[G) + *(H, F) ■ t(H, F)- 1 ■ LF(H). (4) 

Equivalently, for each pair (H, F) considered by the visiting procedure, we 
apply the rule 

|Aut(#UF)| T _ rr . 
X[huf] x [HuF] + | Aut ^| • LF(ff). (5) 

Lemma 3 The total accumulation to x\g] is k ■ LF(G). 

Proof. By Lemma [2] and (jlj), the total accumulation to x\g] from a class 
[(H,F)] satisfying [H U F] = [G] is a(H,F) ■ LF(iJ). Taking the sum over 
all such classes, the claim follows by (jSJ). □ 
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2.3 Implementation details 

Starting with the empty graph of order 14, that is, the unique 0-regular 
graph on 14 vertices, we use the forward accumulation technique for each 
k — 1, 2, . . . , 13 in turn to compute both (i) exactly one graph G from each 
isomorphism class of /c-regular one-factorizable graphs; and (ii) the value 
LF(Cr) for the graphs G. 

The representative graph G in an isomorphism class is the canonical form 
computed by nauty (version 2.2) [11] using the built-in adjtriang vertex 
invariant. The canonical form is stored in 16 bytes of memory as a bit map 
of the upper triangle of the adjacency matrix; the associated accumulator 
variable uses 32 bytes. To enable rapid searching of these 48-byte records, 
we use an open-addressing hash table with 10 9 entries indexed by 4-byte hash 
values of the 16-byte bit maps. 

We use the GNU Multiple Precision Arithmetic Library [5] to carry out 
the accumulator arithmetic. An elementary backtrack algorithm suffices for 
listing the one-factors F G T(H). 

The number of isomorphism classes of fc-regular one-factorizable graphs 
on 14 vertices is, for k = 0, 1, . . . , 13, 

1, 1, 4, 504, 87977, 3459360, 21609293, 21609301, 3459386, 88193, 540, 13, 1, 1. 

The performance bottleneck is at k = 7, where 43,218,594 records need to 
be stored, occupying about 6 GB of memory together with the hash table 
(cf. [15J). In terms of running time, the entire computation (including the 
correctness checks described in what follows) took about 13 days on a Linux 
PC with a 3.66-GHz Intel Xeon CPU and 32 GB of main memory. 

2.4 Correctness checks 

Based on forward accumulation, we have the value LF(G) for every regular 
graph G on 14 vertices. In particular, 

LF(K U ) = 98,758,655,816,833,727,741,338,583,040. (6) 

As a first check, we use Dinitz-Garnick-McKay recursion to verify the 
LF(G) values for each regular G. The same values are obtained both using 
forward accumulation and Dinitz-Garnick-McKay recursion. 

As a second check, we use Meringer's classification program genreg [16] 
to generate all the regular graphs on 14 vertices, and then filter out the 



6 



graphs that do not have a one-factorization. The obtained graphs agree with 
those obtained by forward accumulation. 

As a third check, we use the following observation due to Dinitz, Garnick, 
and McKay. Taking the sum over all isomorphism classes [G] of /c-regular 
graphs on n vertices, we have 

LF(^=("-y'x:^.LF( G ).LF(G). (7, 

Moreover, this holds for every k = 0,1,2, ... ,n — 1. Using data obtained 
from forward accumulation, we evaluate the right-hand side of (|7j) for n = 14 
and k = 0,1, ... ,13. In each case the computed value agrees with (J6j). 
To enable further checks, we display in Table Q] the value 



E 



|Aut(G)l 



14! 

LF(G) 



for each k — 0,1, ... ,13, where the sum is taken over all isomorphism classes 
[G] of /c-regular graphs on 14 vertices. Put otherwise, the tabulated value is 
the number of distinct one-factorizations of /c-regular graphs on a fixed set of 
14 vertices; cf. [HI Table 7], where a slightly different quantity is tabulated, 
however. 



3 Counting nonisomorphic one-factorizations 

Our primary objective in this section is to classify the one-factorizations of 
Ku with nontrivial automorphisms. Once this classification is available, it 
is a simple matter to determine NF(lfi4) based on LF(ifi4) and the clas- 
sification data. Before presenting the classification approach, we discuss a 
representation for one-factorizations of the complete graph in the framework 
of group divisible designs, and narrow down the automorphisms that need to 
be considered to obtain a complete classification. 

3.1 One- factorizations as group divisible designs 

It is well known that a one-factorization of K n can be viewed as a particular 
type of group divisible design (GDD). For our purposes this representation 
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k Distinct one-factorizations 

1 135135 

2 5338373040 

3 78634135419840 

4 461142306338313600 

5 1078882420304271623040 

6 972197327694773750169600 

7 315828427387711768964628480 



8 33491835583595013396417085440 

9 1006698095378044123991615078400 

10 7024525682952576878777802424320 

11 8573318527281503086919968358400 

12 1283862525618838460637401579520 

13 98758655816833727741338583040 



Table 1: Distinct one-factorizations of /c-regular graphs on 14 vertices 

will be convenient as it enables us to immediately prescribe the action of an 
automorphism not only on the vertices but also on the one-factors. 

To develop the GDD representation, let U = {ui,U2, ■ ■ ■ , w„_i} be a set 
with one element for each one-factor, and let V = {i>i,i>2, • • • ,v n } be a set 
(disjoint from U) with one element for each vertex of K n . The elements of 
U U V are called points. We can now let a 3-subset of the form {uk, Vi, Vj} 
carry the information that the edge {vi,Vj} occurs in the one-factor u^. In 
this setting, a one-factorization of K n is a tuple X = (U, V, B), where B is a 
set of 3-subsets of points, called blocks, such that 

(a) for all 1 < k < n — 1 and 1 < i < n, the pair {uk, v occurs in a unique 
block; 

(b) for all 1 < i < j < n, the pair {vi,Vj} occurs in a unique block; and 

(c) for all 1 < k < £ < n — 1, no block contains the pair {uk, ug}. 

A reader acquainted with GDDs immediately observes that X is in fact a 
GDD with group type (n — l) 1 !™, block size k — 3, and index A = 1; however, 
here we find it more convenient to speak of one-factorizations. 
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Two one-factorizations, X = (U,V,B) and X' = (If, V, B'), are isomor- 
phic if there exists a bijection if : U U V — > U' U V such that <p(U) = U', 
(p(V) = V, and f(B) = B'. Such a <p is an isomorphism from X onto X'. An 
isomorphism of ^ onto itself is an automorphism of X. We denote by Aut(X) 
the automorphism group of X. Note that the restriction of an isomorphism 
if to V uniquely determines <p on U. It follows that the standard (graphical) 
and the GDD representations of one-factorizations of K n are equivalent for 
purposes of classification up to isomorphism. 

3.2 Automorphisms of one-factorizations 

Any nontrivial group has a subgroup of prime order, and our classification 
considers all possible such groups that can occur as a group of automorphisms 
of a one-factorization. Ihrig and Petrie [6] carried out an extensive study 
of all possible automorphisms for one-factorizations — in particular, for the 
complete graph Kyi — but we indeed only need those with prime order. 

It is convenient to assume in what follows that the sets U and V are 
arbitrary but fixed. This enables the following two simplifications. First, 
isomorphisms between one-factorizations are permutations of U U V that fix 
U and V setwise. Denote by Y the group of all such permutations of U U V . 
Second, we can identify a one-factorization X with its set of blocks B. 

A group element a G V of prime order p is determined up to conjugation 
in T by the number of fixed points it has in the sets U and V. Denote 
by fu and fy the number of fixed points of a in U and V, respectively. 
Because the cycle decomposition of a consists only of fixed points and p- 
cycles, it is immediate that p must divide both \U\ — fu — n — 1 — fu and 
|V| — fv = n — fv- Not all such types (p, fu, fv) define automorphisms of 
one-factorizations, however. 

We proceed to narrow down the possible types (p, fu, fv)- The following 
lemma is analogous to [3 Lemma 32]; see also [TTl Lemma 4.1]. 

Lemma 4 Let a be an automorphism of a one-factorization X with fu > 
1 and fv > 1. Then, X restricted to the fixed points of a forms a one- 
factorization. In particular, fu = fv ~ 1 an d fv is even. 

Proof. Let x and y be two distinct points fixed by a, at least one of which 
is in V. Such points clearly exist if fu > 1 and fv > 1. Then there is 
exactly one block {x, y, z} of X that contains both x and y. Since a is an 
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automorphism of X, also {x,y,a(z)} is a block, and thus a(z) = z. This 
shows that no block of X intersects the set of fixed points of a in 2 points. 
Consequently, every block of X intersects the set of fixed points in 0, 1, or 
3 points. Disregarding all other blocks except those intersecting in 3 points, 
we obtain a set system X' = (U',V',B'), where U' and V are the sets of 
points fixed by a in U and V, respectively. 

To see that X' is a one-factorization of a complete graph of order n' — fy, 
we first count in two different ways the tuples {uk,Vi,B) such that Uk E U' , 
Vi E V, B E B', and {u k , Vl } C B. We have \U'\ ■ \V'\ = 2\B'\. As all 
2-subsets of V occur in a unique block, we have \B'\ = | V | ( | V | — l)/2. 
Combining the two equalities and noticing that \U'\ = fy and |V'| = fy, we 
get fu = fy — I. Furthermore, it follows that X' meets the requirements (a), 
(b), and (c) in the definition of a one-factorization of a complete graph for 
n' — fy. In particular, fy must then be even. □ 

The following lemma is due to Seah and Stinson [TTj Lemmata 4.2 and 4.3]. 

Lemma 5 Any nonidentity automorphism of a one-factorization satisfies 
fy < n/2. Equality holds only if the automorphism has order 2. 

The following two lemmata rule out certain automorphism types with 
p = 2; the argument in the proof of the second one is similar to that of 
[H Lemma 33]. 

Lemma 6 Let n = 2 (mod 4). If p = 2 and fy — for an automorphism 
of a one- factorization, then fu < n/2. 

Proof. With the objective of bounding fu, consider any one-factor u k E U 
fixed by a. Since the number of edges in the one-factor is n/2 = 1 (mod 2), 
not all orbits of edges in the one-factor can have size 2. Thus, because fy = 0, 
there is at least one block {uk,Vi,Vj} such that Vj) is a 2-cycle of a. The 
number of such 2-cycles in a is n/2, and thereby fu < n/2. □ 

Lemma 7 Let n = 4 (mod 8) or n = 6 (mod 8). If p = 2 and fy — for 

an automorphism of a one- factorization, then fu ^ 1. 

Proof. To reach a contradiction, assume that a is an automorphism of a one- 
factorization with p — 2, f v — 0, and fu = 1. Without loss of generality we 
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Table 2: Automorphism types 

may assume that UUV = { — (n — 1), — (n— 2), . . . ,n— 1} and that = — x 
for all x G U U V. In particular, because a is an automorphism, both U and 
are closed under negation. Because no pair of points occurs in more than 
one block, it follows by a({x, — y, y}) = {a(x),y, —y} that a(x) = x and 
x — 0. In particular, each of the n/2 pairs of the form {— y,y} C V occurs 
in a block of the form {0, —y,y}. All the remaining blocks are moved by a, 
that is, {x, y, z} is a block if and only if {— x, —y, —z} is a block, and each 
such pair of blocks contains either or 4 pairs of points with opposite signs. 
Considering pairs of points with either one point in U and one in V or two 
points in V, there are 2 • (n — 2)/2 • n/2 + {n/2) 2 such pairs of points with 
opposite signs. The fixed blocks account for n/2 occurrences, so the moved 
blocks must account for 3 • n/2 ■ (n/2 — 1) occurrences. For n = 4 (mod 8) 
and n = 6 (mod 8) this number is not divisible by 4, a contradiction. □ 

For n = 14, prime orders p = 2, 3, 5, 7, 11, 13 need to be considered. By 
applying Lemmata 0] to all other automorphisms types (p,fu,fv) except 
those listed in Table [2] are excluded. The last four columns of the table — IF, 
VI, V2, and Seeds — are related to the main search to be discussed in the 
next section. 
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3.3 The classification 



We now present our approach for classifying the one-factorizations that admit 
at least one automorphism of the types in Table [2j For this task we essentially 
rely on the framework established in [7]. 

The classification is based on certain substructures, to be called seeds, 
at least one of which is contained in every one-factorization that we want to 
classify. Due to the requirement of nontrivial symmetry, the precise definition 
of a seed will unfortunately be somewhat technical. A simplified intuition to 
keep in mind is as follows. Consider an arbitrary one-factorization, X, and 
(by some rule) select a set T of points. Let S be the set of blocks in X that 
have nonempty intersection with T. Now, based on the defining properties 
of one-factorizations (and the rule for selecting T), we can anticipate the 
structure of S without knowing all the possible X explicitly, and classify all 
possible S up to isomorphism. These sets S will intuitively be the seeds; the 
added technicality follows because (i) we insist on nontrivial symmetry in 
the form of a prime-order group of automorphisms II < T; and (ii) we must 
make precise the structure of S in relation to IT and T, without reference to 
any containing one-factorization. 

The technical definition of a seed is as follows. Let II < T be a prime-order 
subgroup whose nonidentity elements have one of the types in Table [2j Let 
TC U U V be a set of points such that the type of II and the columns IF, VI, 
and V2 in Table [2] determine the size and composition of T in relation to IT 
as follows. The size of T is determined by the number of columns containing 
either an F ("fixed") or an M ("moved"). The column IF ("one-factor") 
indicates whether T contains a point from U and whether the point is fixed 
or moved by IT. The columns VI ("first vertex") and V2 ("second vertex") 
indicate whether T contains points from V and whether these points are fixed 
or moved by II. (For example, for type p — 3, fu = 1, fv — 2, the set T 
consists of one element of U fixed by II and two elements of V both moved 
by II.) Finally, let S be a union of Il-orbits of 3-subsets of U U V such that, 
referring to the elements of S as blocks, 

(a') for all 1 < k < n — 1 and 1 < % < n, the pair {w^, v{\ occurs in at most 
one block; 

(b') for all 1 < i < j < n, the pair {v j, Vj} occurs in at most one block; 
(c') for all 1 < k < £ < n — 1, no block contains the pair {uk, ue}; 
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(d') for every w fc G T R [/, the point occurs in exactly n/2 blocks; 

(e') for every G T R V, the point occurs in exactly n — 1 blocks; 

(f ') the set T has nonempty intersection with at least one block on every 
Il-orbit on S; and 

(g') the set T occurs in at least one block. 

Each tuple (II, T, <S) meeting these requirements is called a seed. Two 
seeds, (II, T, S) and (II', T", S'), are isomorphic if there exists a 7 G T such 
that 7ll7 _1 = IT, j(T) = T', and 7(5) = 5'. The permutation 7 is an 
isomorphism of (II, T, 5) onto (U',T',S'). An automorphism of a seed is an 
isomorphism of the seed onto itself. We denote by Aut(II, T, S) the automor- 
phism group of a seed. A one-factorization X contains (or extends) a seed 
(n,T,«S) if II < Aut(^) and 5 C ;t. 

We classify the seeds up to isomorphism by enlarging the set T one point 
at a time using the algorithms described in [7] . The number of nonisomorphic 
seeds associated with each automorphism type is displayed in the column 
Seeds in Table EJ 

Because every one-factorization with a nontrivial automorphism group 
contains at least one seed, we can visit every isomorphism class of one- 
factorizations by extending each classified seed in all possible ways. The 
task of finding all one-factorizations that contain a given seed (II, T, S) is 
an instance of the exact cover problem. Put otherwise, we must cover the 
remaining uncovered pairs of points of the form {uk, Vi} and {v j, v j} exactly 
once in all possible ways using Il-orbits of triples of the form {uk, Vi, Vj}; each 
triple covers the pairs that occur in it. For algorithms, we refer to [HI fTUj . 

We reject isomorphs among the visited one-factorizations using the frame- 
work in [7\ , which is an instantiation of the canonical augmentation technique 
developed by McKay [T3]. In essence, we identify a canonical Aut(A')-orbit of 
seeds contained by a visited one-factorization X, and then check whether the 
seed (II, T, S) from which X was extended is in the canonical orbit. If yes, we 
accept X if it is also the (lexicographic) minimum of its Aut(II, T, (S)-orbit. 
Otherwise we reject X . 

In the search we find 581,042,656,543 one-factorizations that extend the 
classified seeds; among these we find 10,305,262,573 nonisomorphic one- 
factorizations with a nontrivial automorphism group. Table |3] displays the 
possible orders i for the automorphism group and the associated number 
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10 


192 


1 


16 


109 


Total 


1132835421602062347 



Table 3: Nonisomorphic one-factorizations of Ku 

Ni of nonisomorphic one-factorizations. (Ni will be determined in Section 
13.41 ) The search was distributed to a network of 180 Linux PCs using the 
batch system autoson [12J, and required in total a little over 5 years of CPU 
time. The prevalent PC model in the network was Dell OptiPlex 745 with a 
2.13-GHz Intel Core 2 Duo 6400 CPU and 2 GB main memory. 

3.4 Applying the Orbit-Stabilizer Theorem 

As a result of the computer searches, we know the value LF(K 14 ) and the 
values % and iV, for each % > 2. We now apply (pQ) in the GDD representation 
over fixed but arbitrary sets U and V. Accordingly, Q is the set of distinct 
one-factorizations over U and V, and T is the group of all permutations of 
UUV that fix U and V setwise. Clearly, |T| = 13! • 14! . Because there are 13! 
ways to label the one-factors in each distinct one-factorization of in the 
standard (graphical) representation, we have = 13! ■ hF(Ku) . Solving 
(JTJ) for Ni and computing J2i>i N, we have that the complete graph K14 has 
exactly 

NF( J ftT 14 ) = 1,132,835,421,602,062,347 
nonisomorphic one-factorizations. 
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3.5 Correctness checks 

We carry out two checks to gain confidence in the correctness of the values 
Ni for i > 2 in Table gfl 

First, both authors independently classified the seeds up to isomorphism, 
with identical results. One classification was conducted using the tools in [7] 
and the other using a basic backtrack search with isomorph rejection based 
on recorded canonical forms. 

Second, for each automorphism type (p, fu, fv), we count in two different 
ways the distinct tuples (X,U,T, S), where X is a one-factorization and 
(II, T, S) is a seed of type (p, fu, fv) contained in X. 

The double count is implemented as follows. Fix an automorphism type 
{p,fu,fv)- For a one-factorization X, denote by seeds(A") the number of 
distinct seeds of the fixed type that X contains. For a seed (II, T, S) of the 
fixed type, denote by ext(II, T, S) the number of distinct one-factorizations 
that extend (II, T, S). Taking the sum over all isomorphism classes of one- 
factorizations on the left-hand side, and over all isomorphism classes of seeds 
of the fixed type on the right-hand side, we have, by the Orbit-Stabilizer 
Theorem, 

in in 

> , ' ' . ■ seeds(^) = V - — ' ' — -. ■ ext(II, T, S). (8) 
jj\Aut(X)\ ^lAut^T^)! 

To evaluate the right-hand side of ([8]), we record the number of extensions 
ext(II,T, S) and |Aut(II, T, S)\ for each classified seed (H,T,S). In partic- 
ular, ext(II, T, S) is simply the number of solutions found in the search for 
exact covers. The right-hand sides obtained for each automorphism type are 
listed in column Count in Table HI 

The left-hand side of OH]) is accumulated for each classified one-factorization, 
X. For each such X, we find all distinct prime-order subgroups II < Aut(X). 
For each such II, we find its type (p,fu,fv), an d accumulate the left-hand 
side of © for this type by |r|/|Aut(Af)| • m(p, fu, f v ), where m(p, fu, fv) is 
the number of distinct seeds of the form (II, T, S) contained in X. Because 
S is uniquely determined by X, II, and T, we can determine m{p, fu, fv) 
by combinatorial arguments based on (p, fu, fv) and Table El For example, 
consider p = 5, fu = 3, fv = 4. From Table [2] we find that T consists of one 
point of U fixed by II and two points of V fixed by IT. By (g') and S C X, 
we have that T is a block of X. There are = 6 possibilities to select a 



15 



V fu fv rn(p,fu,f v ) 



Count 



2 1 

2 3 

2 3 

2 5 

2 5 

2 7 

3 1 
5 3 
7 6 
13 



2 

4 

6 

2 
4 

1 



24 598566905953570569439936512000 

42 10362562621908790673701601280000 

40 11064634216108608459689164800000 

70 109764651070947200382428774400000 

48 314439007330643189170176000000 
98 65958792604977492862823301120000 

1 814728009186504568995840000 

6 1840950333789938122752000 

49 2850020421206999040000 
13 6016709778103664640000 



Table 4: Double counting check 



pair of points in V fixed by II, each of which occurs in a unique block of X . 
Thus, m(5, 3, 4) = 6. The other values m(p, fu, f v ) are displayed in Table HI 
For each type (p, fu, fv), we find that the computed left-hand and right- 
hand sides of ([HD agree. This, together with the observation that the left-hand 
side of ([8]) depends on each of the computed values % and Ni for i > 2, gives 
us confidence that Table [3] is correct. 
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